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, Abstract. A high quality, medium-resolution K-band spectrum has been obtained of the ultracompact Hll re- 

' gion G29.96— 0.02 with the Very Large Telescope (VLT). The slit was positioned along the symmetry axis of 

^Nj ' the cometary shaped nebula. Besides the spectrum of the embedded ionizing O star, the long-slit observation 

, revealed the rich emission-line spectrum produced by the ionized nebula with sub-arcsec spatial resolution. The 

' nebular spectrum includes Br7, several helium emission lines and a molecular hydrogen line. A detailed analysis is 

^ I presented of the variation in strength, velocity and width of the nebular emission lines along the slit. The results 

are consistent with previous observations, but the much better spatial resolution allows a critical evaluation of 
j ' ' models explaining the cometary shape of the nebula. Our observations support neither the wind bow shock model 

^ , nor the champagne flow model. 

^— ■») ' The measured line ratios of the nebular hydrogen and helium lines are compared to predictions from case B 

^■f^ ' recombination-line theory. The results indicate an electron temperature between 6400 and 7500 K, in good agree- 

, ment with other determinations and the Galactocentric distance of 4.6 kpc. The He'''/H+ ratio is practically 

constant over the slit; we argue that He is singly ionized throughout the nebula. We review the various observa- 
tional constraints on the effective temperature of the ionizing star and show that these are in agreement with its 
K-band spectral type of 05-06 V. 
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1. Introduction The frequently seen cometary structures are suggestive of 

j; . bow shocks. Wood & Churchwell (1989) proposed that 

. H II regions trace the formation sites of massive stars. The ^^^^^ ^ow shocks might resuh from the supersonic motion 

ionizing radiation of the embedded, young OB stars is ab- ^ ^ kms"!) of the embedded O stars relative to 

sorbed by gas and dust in their near surroundings and ^.^e ambient molecular cloud. The stand-off distance of the 

re-emitted at infrared and radio wavelengths. Their pro- ^^^^^^ ^^^^ observed "size" of the UCHii regions, 

genitors, the ultracompact Hii (UCHii) regions, represent determined by the density of the ambient medium, the 

the earliest recognizable phase of massive-star formation, velocity of the O star, and the momentum of its stellar 

and are the most luminous sources at 100 /im, observ- ^-^^^ ^hese parameters are not expected to vary rapidly 

able throughout the Galaxy (Churchwell, 1990; Garay & ^j^h time, thus delaying the apparent expansion rate of 

Lizano, 1999). The small size and large number of UCHii ^-^^ uCHil region (e.g. van Buren et al., 1990; Mac Low 

regions indicate that the lifetime of this phase is signif- ^gg^. ^^^^^ ^ Low, 1992). 
icantly longer than the ~ 10"* year predicted by models 

describing the expansion of young Hii regions. G29.96— 0.02, also denominated IRAS 18434—0242 

A possible solution for this "confinement problem" is (hereafter G29.96), is a well-studied, cometary-shaped 

suggested by the observed morphologies of UCHli regions. UCHii region at a Galactocentric distance of 4.6 kpc (e.g. 

Pratap et al., 1999). Arguments based on extinction and 



Send offprint requests to: N.L. Martin-Hernandez, email: leti- tl^e spectral type of the ionizing star favour the near helio- 

cia.martin@obs.unige.ch centric distance of 6 kpc instead of the far distance of 9 kpc 

* Present address: Geneva Observatory 1290 Sauverny (Pratap et al., 1999). The cometary structure of G29.96 

Switzerland was first observed by Wood & Churchwell (1989) at radio 
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wavelengths. Follow-up radio observations were performed 
by Wood & Churchwell (1991), Afflerbach et al. (1994), 
Cesaroni ct al. (1994), Fey ct al. (1995) and Kim & Koo 
(2001). The cometary shape is characterized by a sharp, 
bright, arc-like leading edge and a low-surface brightness 
tail of emission which trails off away from the edge. Wood 
& Churchwell (1989) found the bow shock hypothesis con- 
sistent with their radio recombination line observations 
of G29.96. Alternatively, a champagne-flow model, which 
describes the spherical expansion of an H ii region into a 
molecular medium with a steep density gradient (Yorke, 
1986), has been proposed to explain the morphology and 
kinematical structure of G29.96 (Lumsden & Hoare, 1996, 
1999). 

G29.96 is embedded in a molecular cloud that pro- 
duces ~ 15 magnitudes of visual extinction (see Martin- 
Hernandez et al., 2002a). This molecular cloud has been 
extensively studied at low spatial resolution (Churchwell 
et al., 1990, 1992). Strong emission from different dense 
gas tracers indicates the presence of hot, dense gas. Just 
in front of the edge of the Hll region, at approximately 
2" west of the radio continuum peak, a hot molecular core 
was discovered in interferometric observations (Cesaroni 
et al., 1994, 1998; Maxia ct al., 2001). Recently, this hot 
core has been found at mid-infrared wavelengths too (De 
Buizer et al., 2002). The observed offset between the hot 
core and the Hii region and the detection of masers to- 
wards the core (Pratap et al., 1994; Hofner & Churchwell, 
1996; Walsh et al., 1998) suggest that the hot core is not 
heated by the ionizing star of G29.96, but by another 
young, deeply embedded massive star. 

G29.96 seems to be associated with a small cluster of 
stars (Fey et al., 1995; Lumsden & Hoare, 1996; Watson 
et al., 1997; Pratap ct al., 1999). The main ionizing star is 
located in the focus of the bright arc. Watson & Hanson 
(1997) obtained a K-band spectrum of the embedded ion- 
izing star, the first ionizing star of an UCHii region ever 
detected. Based upon the presence and strength of pho- 
tosphcric lines of C iv and N ill, they deduced a main se- 
quence spectral type between 05 and 08. 

G29.96 provides an excellent laboratory to test mod- 
els describing the ionization structure of an UCHii region. 
The ionization structure can be traced by the nebular hy- 
drogen and helium recombination lines. The comparison 
between observations and model provides an important 
constraint on the stellar parameters of the ionizing star(s), 
yielding important information on the properties of newly 
formed massive stars (e.g. Vanzi et al., 1996; Lumsden & 
Puxley, 1996; Lumsden et al., 2001a; Hanson et al., 2002). 
Due to the large visual extinction towards H II regions like 
G29.96, such a study has only recently become possible 
thanks to the technological development of near-IR spec- 
trometers. 

In this paper we present new, high quality, long-slit 
K-band spectra of G29.96 (2.07-2.19 fim). Besides Br7, 
a large number of He i lines and a H2 line are covered. 
The study of the relative strengths of the Hi and Hei 
lines as a function of position along the slit allows differ- 



ent lines of investigation. First, it is possible to study the 
impact of large line opacities on the Hei lines. Second, the 
conditions present in the ionized gas can be determined. 
Third, the results constrain the effective temperature of 
the ionizing star, which can be compared with the spec- 
tral type implied by its K-band spectrum. Finally, the 
kinematical structure of the ionized gas can be studied at 
higher spatial resolution than in previous studies (Wood 
& Churchwell, 1991; Afflerbach et al., 1994; Lumsden & 
Hoare, 1996, 1999) and compared with the molecular gas 
motions. This comparison can be used to test the different 
models (bow shock, champagne flow, etc.) that have been 
proposed to explain the cometary morphology of G29.96. 

This paper is structured as follows. Section 2 de- 
scribes the observations and the data reduction. Section 3 
presents the resulting spectra and the variations and kine- 
matics of the lines across G29.96. Section 4 describes the 
H I and He i recombination theory and shows the differ- 
ent applications of these lines. Section 5 reviews the vari- 
ous observational constraints on the effective temperature 
of the ionizing star. Section 6 discusses the dynamics of 
G29.96. Finally, Sect. 7 summarizes the results of the pa- 
per. 

2. Observations 

Long-slit (120") K-band spectra were obtained of the 
UCHii region G29.96 with ISAAC mounted on Antu 
(UTl) of ESO's Very Large Telescope (VLT), Paranal, 
Chile. The observations were carried out on 19 March 
2000. A slit-width of O'.'S was used, resulting in a spec- 
tral resolving power R~ 8000 (=c/Av). The slit was po- 
sitioned along the symmetry axis of the object (116.2 de- 
grees with respect to the north-south axis). Figure 1 shows 
the location of the slit in the (acquisition) image of G29.96 
obtained through a narrow-band H2 filter (central wave- 
length 2.13 /im). 

The observing conditions at Paranal Observatory were 
excellent (humidity less than 10% and seeing of 0'.'6). In 
order to correct for the sky backgromid, the object was 
"nodded'' between two positions on the slit (A and B) 
such that the background emission registered at position 
B (when the source is at position A) is subtracted from 
the source plus sky background observation at position B 
in the next frame, and vice versa. This strategy has the 
advantage that observing time is very efficiently used, but 
it reduces the effective length of the slit by about 50%; for 
G29.96 this is not a problem. 

The electrical ghosts and bias were removed from the 
frames before flatfielding. Subsequently, the sky back- 
ground emission was subtracted following the procedure 
outlined above. Wavelength calibration was performed us- 
ing the telluric OH lines. The accuracy of the wavelength 
calibration is 3 kms~^. To correct for absorption lines 
due to the Earth's atmosphere, the target spectra were 
divided by the spectrum of a telluric standard star. For 
this purpose, an A-type star was observed under identical 
sky conditions, which provides a continuous spectrum only 
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Fig. 1. Acquisition image of the UCHll region 
G29.96-0.02 obtained through a narrow-band H2 
1-0 S(l) filter (central wavelength 2.13 /xm). 

Coordinates arc relative to the position of the 
ionizing star; a(2000) = 18'' 46^4^075, 5(2000) = 
-2°39'21'.'42. Contours are 1, 2.5, 5, 10, 25, 50, 75 
and 90% the peak flux. The cometary structure of 
the Hii region is clearly revealed. The position of 
the slit is marked by a dashed lino. The ionizing star 
and the bright arc are indicated. At a distance of 
6 kpc, 1 arcsec corresponds to 0.029 pc. 
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containing telluric absorption lines. The hydrogen Br7 line 
in the A- star spectrum was removed by a model fit. The 
telluric standard was also used to correct for the through- 
put of t(^l(^scopc> and instrument. For a more detailed de- 
scription of the data reduction procedures we refer to Bik 
et al. (in preparation). 

Spectra were extracted along the slit (1 pixel = 0'.'147) 
in the spatial direction. The peak positions, equivalent 
widths and fluxes of the lines were measured by fltting a 
Voigt profile. A line is defined as being detected if its peak 
intensity exceeds the rms noise level of the local continuum 
by at least a factor three. The central wavelength of the 
line was measured with respect to the local standard of 
rest, i.e. corrected for the projected velocity of the Earth 
with respect to the Sun at the time of observation and 
for the solar motion in the direction of G29.96 (fLSR = 
44 kms~^). 

3. Results 

Figure 2 shows the nebular spectrum of G29.96 extracted 
from a O'.'S x 0'.'882 region around the peak of the emission. 
The strongest line corresponds to Br7. The detected He I 
emission lines are the transitions of 4^S-3^P at 2.1128 /zm, 
4IS-3IP at 2.1137 /im, 73S-43p at 2.1501 /xm, 7^Y-A^V> 
at 2.1615 Atm, 7^Y-A^V) at 2.1624 urn, the I'^^^G-A^^^Y 
blend at 2.1649 /um and the transitions of 7^P-4^D and 
7ID-4IP at 2.1824 and 2.1848 /xm, respectively. From 
these Hex lines, significant emission along the slit was 
found for the transitions of 43S-3^P at 2.1128 /xm, A^^- 
3ip at 2.1137 /xm, 7^Y-A^T) at 2.1615 /im and the l^^^G- 
4i'3F blend at 2.1649 /tm. The 1-0 S(l) transition of H2 at 
2.1218 /xm and the forbidden line of [Feiii] at 2.1457 /xm 



arc also evident. Table 1 lists the observed lines with their 
identifications and their fiuxcs relative to that of Br7 at 
the peak position. An unidentified feature is present at 
2.1022 /im. This feature has also been detected in near- 
infrared spectra of compact planetary nebula (Lumsden 
et al., 2001b). 

We could measure the nebular line emission along 
~ 20" of the slit, centered on the position of the ioniz- 
ing star. Figure 3 illustrates the spatial variation of the 
strongest lines in this region. Moving from right to left 
along the slit (see Fig. 1), the Br7 flux (Fig. 3a) slowly 
increases until the position of the star. After the star, the 
Br7 flux rapidly reaches a maximum and subsequently de- 
creases. This maximum is at 2'.'27 ahead of the star and 
coincides with the arc observed at radio wavelengths (e.g. 
VVood & Churchwell, 1989). The peak corresponding to 
the arc is fitted by a Gaussian profile with a full width at 
half maximum (FWHM) of 1'.'14. The global structure of 
the tail to the right of the star can be fitted by a semi- 
Gaussian centered at the position of the star and with a 
FWHM of 7'.'1 and a height 0.21 times the peak fiux. The 
following panels (b f) in Fig. 3 compare the flux distri- 
bution of the strongest Hel lines and the H2 2.1218 /tm 
line to that of Br7. The Hex and Br7 distributions trace 
each other very well, indicating that the zones of H+ and 
He+ are of roughly the same size. This indicates that the 
central ionizing star must be hot enough to produce ion- 
ized helium throughout the Hix region (see Sect. 4.2 for a 
detailed discussion). The variation in the H2 line strength 
resembles the ionized gas profile well, although it displays 
some significant deviations. It peaks 0'.'93 ahead of the Br7 
maximum and exhibits a different structure in the tail of 
G29.96. The H2 peak corresponding to the arc can be fit- 
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2.14 

Wavelength (/im) 

Fig. 2. Nebular spectrum of G29.96-0.02 extracted from a O'.'S x 0'.'882 region around tlie peak of ionized emission. The lines 
observed, their identification and their relative flux with respect to Br7 are listed in Table 1. Note the logarithmic scale of the 
flux scale. 



Table 1. Emission lines in the nebular K-band spectra of G29.96— 0.02. 



At (urn) 


Identification 


Peak flux* 


At (Aim) 


Identification 


Peak flux* 


2.10213^ 


? 


0.004 


2.16147 


He 73F-43D 


0.024 


2.11282 


He 433-33? 


0.041 


2.16240 


He 7IF-4ID 


0.007 


2.11371 


He 4IS-3IP 


0.011 


2.16495 


He 7i'3G-4i'3F 


0.046 


2.12183 


H2 1-0 S(l) 


0.005 


2.16612 


H 7-4 (Br7) 


1.00 


2.14570 


[Feiii] 3G3-3H4 


0.009 


2.18245 


He 7^P-A^D 


0.005 


2.15007 


Ho 7'^S 4-'^P 


0.002 


2.18483 


He 7ID 4ip 


0.001 



(t) Vacuum wavelength. The vaccum rest wavelength for the H2 line is taken from Bragg et al. (1982) and those of the Hei lines 
from Benjamin et al. (1999). (*) Flux with respect to Br7 at the position of the arc peak. (^) Observed line center wavelength. 



ted by a Gaussian profile with a FWHM of 0'.'87. Since 
the spatial profiles of the Hi, He i and H2 lines have been 
obtained simultaneously in one same setting, the relative 
separation of the Br7 and H2 peaks is reliable. 

The velocity structure of the ionized (as traced by, for 
instance, Br7) and molecular gas (as traced by H2) along 
the symmetry axis of G29.96 is shown in Fig. 4. Figure 4 a 
plots the centroid velocity of the Br7 line along the central 
20" of the slit. There is a large variation in radial velocity 
(~ 38 km s~^) across the object. The gas in front of the 
head of the arc has a velocity which increases from 67 (po- 
sition —7") to 87 kms~^ (position —3"). At the arc, be- 
tween positions —3" and —1", the velocity increases more 
rapidly, reaching a value of 99 kms~^ at the position of 
the line intensity peak. It continues to increase towards 
the tail where it reaches 105 kms~^ at approximately 3" 
behind the star. Afterwards, the velocity drops to a value 
of approximately 94 kms"""^ at ~ 5" behind the star. 

Figure 4 c shows the behaviour of the Br7 FWHM. 
As can be seen, there is a very rapid increase of the 
linewidth ahead of the line intensity peak, which gets up 
to ^ 62 kms~^ wide. The Br7 line is ~ 42 — 46 kms~^ 
wide behind the arc. 



This behaviour of the Br7 line velocity across the sym- 
metry axis of G29.96 agrees well with the radio recom- 
bination line studies of Wood & Churchwell (1991) and 
Afllerbach et al. (1994). Low spatial resolution observa- 
tions of Br7 have been presented by Lumsden & Hoare 
(1996, 1999). In contrast to their results, we see a larger 
variation in the peak velocity of the Br7 line along the 
symmetry axis. In particular, the appreciable increase in 
velocity by more than 30 kms~^ just ahead of the arc is 
not apparent from their study. This is likely due to the 
poorer spatial resolution of this earlier study combined 
with the steep increase in peak intensity towards the arc, 
which gets a high weight in their spatial averaging proce- 
dure. 

The velocity structure observed for Br7 is compared 
to that of the molecular gas traced by the H2 1-0 S(l) 
line (see Figs. 4 b and 4d). There is little variation in the 
peak veloc'ity of the H2 line, which shows an almost con- 
stant value of 100 kms~^ across G29.96 (calculated be- 
tween positions —3" and -1-5", where the H2 line emission 
is significant, i.e. above 3a). The average velocity of 100 
kms~^ is consistent with the values (around 98 kms~^) 
obtained from single dish observations of rotational tran- 
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Fig. 3. Spatial variation of a) Br7 , b— e) the strongest Hex recombination lines and f) H2 1-0 S(l) across the symmetry axis of 
G29.96 (from left to right in Fig. 1). The normalized He I and H2 intensity distributions are compared to the Br7 profile (dashed 
lino). The origin corresponds to the position of the ionizing star. At a distance of 6 kpc, 1" corresponds to 0.029 pc. While the 
He I and Hi distributions trace each other very well, the H2 line peaks 0'.'93 ahead of Br7 and shows a different structure in 
the tail. We note that it was not possible to measure line fluxes close to the position of the ionizing star, which explains the 
absence of points in between —1" and 0". 



sitions of CS (e.g. Churchwell et al., 1992; Bronfman et al., 
1996; Olmi & Cesaroni, 1999) and NH3 (e.g. Churchwell 
et al., 1990; Cesaroni et al., 1992). VLA measurements of 
ammonia by Cesaroni et al. (1994) resulted in a value of 
98.7 kms~^ towards the hot core detected just in front 
of the arc of G29.96. Thus, the H2 gas is practically at 
rest with respect to the ambient medium. The H2 FWHM 
shows a practically constant value of 33 kms^^ in the in- 
terval [0", +7"]. This FWHM is of the order of the spectral 
resolution, implying that the H2 line is unresolved between 
these positions. Similarly to the case of the ionized g as, a 
significant increase of the FWHM up to ^ 50 kms~^ is 
observed in front of the star. 

The K-band spectral observations of G29.96 by 
Lumsden & Hoare (1999), although of much lower spatial 
resolution than ours, have about twice our spectral reso- 
lution (~ 18 kms~^). Their observations of the 1 — S(l) 



H2 line show, as well, a very rapid increase in linewidth 
(from ~ 20 to 28 kms~^ in their case) ahead of the arc. 
They were also able to discern small variations of the H2 
centroid velocity in the same direction as those of the ion- 
ized gas. This behaviour is, however, not evident from our 
data. 

Comparing the Br7 and H2 velocity structure, we note 
that around the position of the ionizing star (between posi- 
tions —2" and +3"), the ionized gas is systematically red- 
shifted relative to the molecular gas by a few (2-5) km s^^ , 
whereas in the tail and ahead of the arc it becomes pro- 
gressively more blueshifted. The increasing blueshift of the 
ionized gas in front of the arc coincides with the signifi- 
cant increase in the Br7 and H2 FWHM. This pattern has 
to be due to acceleration of the ionized gas along the line 
of sight. 
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Fig. 4. Spatial variation of the velocity centroid (top panel) and line FWHM (bottom panel) of Br7 and H2 along the symmetry 
axis of G29.96 (from left to right in Fig. 1). Typical error bars are indicated at different positions along the slit. In panels a) 
and c), the intensity distribution of Br7 is plotted (dashed line) for comparison. The dotted lines in panels b) and d) represent 
the average H2 centroid velocity and FWHM calculated between positions —3" and +5", and 0" and +5", respectively. The 
zero point along the slit corresponds to the position of the ionizing star. At a distance of 6 kpc, 1" corresponds to 0.029 pc. 



4. The recombination lines 

The nebular spectrum of G29.96 is dominated by recom- 
bination lines of hydrogen and helium. For most of the 
Hi emission lines observed in nebulae, there are no radia- 
tive transfer problems; hence, the Case B approximation 
(Hummer & Storey, 1987), which assumes that the Lyman 
transitions are optically thick, reproduces well, apart from 
dust extinction, the observed nebular spectra. The recom- 
bination radiation of He l singlets is very similar to that 
of Hi, and can also be approximated by the Case B treat- 
ment. In this case, by considering Case B it is assumed 
that the nebula has a large optical depth in transitions 
arising from the l^S level. However, the recombination 
cascade of the Hci triplets can be modified by the fact 
that the 2^S level (see Fig. 5) is highly metastable. All 
recaptures to triplets eventually cascade down to the 2^S 
level. Depopulation of this level can only occur through 
photoionization, through collisional transitions to 2^S and 
2ip, or through the forbidden 238-1^8 transition. The 



rates of these processes are small and hence, the popu- 
lation of the 2'^S level is large; in turn, the optical depth 
of the lower n^P-2^8 transitions is substantial. Moreover, 
the collisional transitions from to upper levels become 
significant and the emission from these levels is enhanced 
with respect to the predictions of a pure recombination 
model. Because of this interplay of collisional and radia- 
tive transfer effects, the triplet lines become sensitive to 
the density and optical depth. 

Recently, Benjamin et al. (2002) have considered the 
effect of the optical depth of the 2*^8 level on the recombi- 
nation spectrum of a spherically symmetric nebula. They 
used a model atom with individual levels up to n = 20 and 
parameterized their calculations in terms of the line center 
optical depth of the 33p-238 3890 A line. A FORTRAN 
program to calculate emissivities of lines arising from lev- 
els with n < 10 over a large range of nebular conditions is 
available from Benjamin et al.^ 

^ http://wisp.physics.wisc.edu/~benjamin 
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Fig. 5. Partial Grotrian energy-level diagram of the He I singlet and triplet transitions showing the infrared lines observed in 

the K-band spectrum of G29.96 (solid arrows), other important infrared lines (dotted arrows) and the strongest optical linos 
(dashed arrows). Note that the energy levels are not drawn to scale. The wavelengths of the infrared lines are indicated in /xm; 
those of the optical lines in A. The wavelengths of the lines used in the present analysis are marked by a box. 



1.8 



10 


1 

Hel 2.1128 _ 


50 




30. 




15, 

lol 

5 1 




2 i 
J 




II 



5000 8000 10000 15000 

Te (K) 

Fig. 6. Optical depth correction factor, /, for the 2.1128 /xm 
Hel line as a function of the electron temperature for rasgo 
between and 100 and an electron density rie = cm~^. 

We use this program to estimate the Case B line fluxes 
of the Hel lines observed in the spectrum of G29.96 and 
their possible enhancement due to coUisional excitation 
from the 2^S level, or as a result of significant opacity in 
the n'^P scries. 

Although Benjamin et al. (2002) include the coUisional 
coupling between the triplet and singlet levels, they find 
that only the triplet lines arising from low levels show a no- 
ticeable enhancement in strength with respect to the zero 
optical depth case. From the triplet lines observed in the 
spectrum of G29.96, wc find that the effect of the finite op- 
tical depth of the 2*^8 level is only evident in the transition 
of 4^S-3^P at 2.1128 /xm. Adopting the nomenclature by 



Benjamin ct al. (2002), wc define the optical depth correc- 
tion factor /lino("-e,re,T389o)= j (^o ,T3890 ) /j (^e ,0) , 

i.e. the ratio of the line emissivity for optical depth raggo 
to the emissivity for zero optical depth. Figure 6 presents 
the variation of the optical depth correction factor for the 
2.1128 /im line as a function of raggo and Te for the case 
of a nebula with rio^lO^ cm^'^ 

The optical depth of the 3^? 2'''S 3890 A fine in 
G29.96 can be roughly estimated. We have that T3890 = 
n(2'^S)(T389o.R, where n(2^S) is the density of atoms in the 
metastable state, (T3890 is the cross section at the line cen- 
ter of the 3890 A line and R is the radius of the nebula. 
The relative population of the 2^S level depends on the 
local electron density and temperature, and can be es- 
timated using equation 5 in Kingdon & Ferland (1995). 
Assuming an electron density of 2x10^ cm~^ and an elec- 
tron temperature of 6500 K (Afflerbach et al., 1994), 
we find that n(23S)/n(He+) - 8 x 10"*^ in the tafl of 
G29.96. Equation 1 in Benjamin et al. (2002) gives that 
a-3890 ~ 7 X 10"" cm2 for Te= 6500 K and a turbulent 
velocity of km s~^. Assuming a helium abundance of 
10% by number and a nebular radius of 7" (sec Sect. 3), 
which at a distance of 6 kpc corresponds to ~ 0.20 pc, we 
obtain T3890 ^ 100. Considering a turbulent velocity of 
20 kms^^ and a filling factor of 0.1, we obtain T3890 ~ 80. 
Hence, the 2.1128 jjum line is expected to be enhanced (see 
Fig. 6) with respect to the pure recombination model. 

It is illustrative at this point to use the relative 
strengths of the Hel recombination lines to (1) determine 
how closely these lines follow the recombination theory 
and (2) study their dependence on the local nebular con- 
ditions. The ratio between any two Hei lines is insensitive 
to the electron density, but can depend on the electron 
temperature. 



8 



N.L. Martm-Hernandez et al.: A VLT spectroscopic study of the ultracompact Hii region G29.96— 0.02 




2468 2468 
Hel 2.1649/HeI 2.1137 Hel 2.1649/HeI 2.1137 

Fig. 7. a) Plot of the line ratio 2.1615-7-^F 4''D/2. 1137-4^8 3^P against the ratio 2.1649-7^ -^G 4^ '^F/2. 1137-4^8 3^P. Note 
that these hne ratios should not depend on radiative transfer effects, b) 2.1649-7^'-^G-4^ -^F/2.1128-4'^S -3'^P against 2.1649- 
7i,3Q 4i.3p/2.1138-4^S-3^P. In both figures, the solid line represents the Case B predictions for Tc = 5, 7.5, 10 and 15 xlO^ K 
(indicated by solid diamonds) and rie = 10* cm~^. The data indicate a temperature between 5000 and 7500 K. The non-solid 
lines in b) represent the calculations by Benjamin et al. (2002) for non-zero optical depth of the 2^S metastable level. Only 
points located in the interval [—5", +8"] are considered. The data points from the tail of G29.96 (> 0") are plotted as open 
circles; those from the arc (< 0") are plotted as solid circles. The open and filled circles have very similar distributions. 



4.1. The electron temperature 

Figure 7 a shows the relation between the Hei line ra- 
tios 2.I6I5-73F-43D/2.II37-4IS-3IP and 2.1649-7i'3G- 
4i'3F/2.1137-4iS-3ip. These ratios show the largest de- 
pendence on Te- The calculations by Benjamin et al. 
(2002) show that these ratios are expected to be inde- 
pendent of radiative transfer problems arising from the 
metastable 2'^S level. The solid line represents the theo- 
retical predictions by Benjamin et al. (2002) for Te = 5, 
7.5, 10, 15 xlO^ K and no = 10^ cm'^. The model re- 
produces the trend in the data reasonably well, although 
there is a slight mismatch between data and theory. The 
data points are, on average, 20% below the model predic- 
tions. The origin of this discrepancy is unclear. It is clear 
from Fig. 7 a that the electron temperature of the ionized 
gas is in between 5000 and 7500 K, and that both the 
points from the tail of G29.96 (> 0") and those from the 
arc (< 0") are distributed equally. 

Figure 7 b shows the relation between the He i lino 
ratios 2.1649-7i'3G-4i'3F/2.1128-43S-33p and 2.1649- 
7i,3Q_4i,3p/2.ii37_4iS-3ip. The data show a significant 
spread in these line ratios, largely exceeding the quoted 
uncertainties. We conclude that the 2.1649/2.1128 /xm 
line ratio is sensitive to the optical depth and that the 
2.1128 /Ltm line is enhanced with respect to the pure 
recombination case (solid line). A comparison with the 
model calculations (Benjamin et al., 2002) indicates opti- 



cal depths, T3890, in between and 100, as expected (see 
above). 

A detailed electron temperature profile along the axis 
of G29.96 can be obtained by using the ratio of any two 
of the He I recombination lines. Figure 7 a shows that the 
2.1649/2.1137 /xm and 2.1615/2.1137 /xm line ratios are 
reasonably well reproduced by the Case B approximation 
and can be used to probe variations in Tc. Hence, by com- 
paring the observed line ratios at any position with the 
Case B predictions we can measure the variation of 
along the slit. 

Figure 8 a shows the Te profile along the symmetry 
axis of G29.96 as derived from the 2.1649/2.1137 /xm line 
ratio. We see that the electron temperature is, within the 
uncertainties, constant along the slit with a mean value of 
6400 ± 150 K. When we use the 2.1615/2.1137 /zm line ra- 
tio, we also obtain a similarly constant To structure along 
the slit, but with a higher average value of 7550 ± 150 K. 
The difference between these two determinations comes 
from the discrepancy between data and theory shown in 
Fig. 7 a. 

Our measurement of Te is in good agreement with pre- 
vious determinations. The non-LTE analysis of radio re- 
combination line observations by AfHerbach et al. (1994) 
resulted in electron temperatures in between 6200 and 
8600 K. They estimate an average Te of ~ 6500 K. The 
electron temperatures derived by Wood & Churchwell 
(1991) from H76q; recombination line observations to- 
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Fig. 8. a) Electron temperature structure along the symmetry 
axis of G29.96 (from left to right in Fig. 1) as derived from the 
He I line ratio 2.1649/2.1138 /um. The dashed line represents 
the average Tc of 6400 K. b) Hc+/H+ structure. The dashed 
Une represents the average value of 0.087. Only points located 
in the interval [—5", +8"] are considered. 



wards G29.96 are low, ranging from ~ 2500 K in the tail 
to 4200 K in the arc. However, for this line, deviations 
from LTE can lead to temperature estimates lower than 
the true nebular temperature. Their best estimate of the 
electron temperature using non-LTE models is ~ 5000 K 
throughout the nebula. Finally, Watson et al. (1997) ob- 
tained temperatures in between 5000 and 7000 K over the 
nebula using 2, 6 and 21 cm continuum data. 

Our estimate of To is also consistent with the 
Galactocentric distance of G29.96. The electron tempera- 
ture in Hll regions has been shown to increase with the 
distance to the Galactic Center as a consequence of the 
decrease in metal content. This To gradient is approxi- 
mately given by Te=5000-h5000xi?Gai/15 (Shaver et al., 
1983; Afflerbach et al., 1996; Deharveng et al., 2000). At 
the Galactocentric distance of G29.96 (4.6 kpc), the ex- 
pected Te is 6500 K, consistent with our value. This elec- 
tron temperature corresponds to an enhancement of the 
metal content by a factor of ~ 2 with respect to that in 
the solar neighbourhood. This metallicity is in agreement 
with the analysis based on the infrared fine-structure lines 
observed towards G29.96 (see Martin- Hernandez et al., 
2002a). 



4.2. He+/I-I+ 

The He+/H+ structure along the symmetry axis of G29.96 
can be determined from the ratio of a singlet He i line to 
Br7 once a specific model for the To structure of the H ii 
region has been assumed (such a ratio is insensitive to 
the electron density). The only singlet line we were able 
to measure along the slit is the transition of 4^S-3^P at 
2.1137 fxm. 

Figure 9 shows the relation between this ratio and 
the Te-indicator Hei 2.1649/2.1137 fim. There is an 
(anti) correlation between these two line ratios because of 
the different Te dependence for the Hei single/triplet and 
the H I lines. The observed data points are well matched by 
He+/H+ abundances in between 0.08 and 0.10 and elec- 
tron temperatures in between 5000 and 7500 K. A pre- 
ponderance of low He"'"/H+ abundances in the tail region 
is evident from this figure. This may reflect a He+ zone 
somewhat smaller than the H+ zone. Such a difference 
would be spatially more prominent in the low density tail. 

The detailed He"'"/H+ structure along the symmetry 
axis of G29.96, calculated for Te = 7000 K, is plotted 
in Fig. 8 b. We find that He+/H+ is practically constant 
across the slit. We obtain an average He"'"/H+ abundance 
of 0.087±0.001. Considering a total helium abundance by 
number of ~ 0.1, this implies that He is practically singly 
ionized throughout the nebula. 

Kim & Koo (2001) estimated the value of He+/H+ at 
different positions in G29.96 from the ratio of the ratio 
of the radio recombination lines He76a and H76q:. They 
derived values in between 0.063 and 0.080, slightly lower 
than our estimate of 0.087. Their lower estimate could be 
due, however, to non-LTE effects on the radio recombina- 
tion lines that may not have been taken into account, or 
to beam smearing. 

We have already mentioned in Sect. 3 that the central 
ionizing star in G29.96 must be hot enough to produce 
ionized helium throughout the H II region. In order to in- 
vestigate the dependence of the relative extent of the He"*" 
region compared to that of the H+ region, Rue+fRa+j 
on the effective temperature of the star, we have calcu- 
lated a set of nebular models with the photoionization 
code CLOUDY version 94.00 (Ferland 2000^) using MICE, 
the IDL interface for CLOUDY created by H. Spoon^. We 
compute nebular models for a static, spherically symmet- 
ric, homogeneous gas distribution with one ionizing star in 
the center. An inner cavity with a radius equal to 10^^ cm 
is set. Two different sets of stellar atmosphere models are 
taken to describe the spectral energy distribution (SED): 
the CoStar models by Schaerer & de Koter (1997) and 
the WM-Basic models by Pauldrach et al. (2001). We use 
stellar models for main sequence (dwarf) stars. The num- 
ber of hydrogen ionizing photons emitted by the central 
source is fixed to that of the SED used, the metallicity of 
the nebula is set to twice solar and the density equal to 



^ see also http://thundor.pa.uky.cdu/cloudy 
^ http://www.astro.rug.nl/~spoon/mice.html 
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Fig. 9. The hnc ratio Hel 2.1649/Br7 is plotted against the 
ratio Hel 2.1649/2.1137. The lines represent the theoretical ra- 
tios for 3 different He+/H+ abundances (0.08, 0.10 and 0.12) 
and an electron density of 10'* cm~^. The values correspond- 
ing to Te = 5, 7.5 and 10 xlO^ K are indicated by solid dia- 
monds. As in previous figures, only points located in the in- 
terval [—5", -1-8"] are considered, the data points in the tail of 
G29.96 are plotted as open circles and those from the arc as 
solid circles. 
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Fig. 10. The relative extent of the He"*" region to that of the 

H+ region, 7?^^+ /7?jj+ , is plotted against Tcff, the effective 
temperature of the ionizing star. For Tes > 40 000 K, the He"*" 
and H"*" zones are coincident. 



10^ cm"^, appropriate for G29.96 (AfHerbach et al., 1994; 
Wood & Churchwell, 1991). 

Figure 10 shows i?Hc+/^H+ as a function of the ef- 
fective temperature of the star, Tes- It is seen that for 
TeS > 40 000 K, the He+ and H+ zones are coincident, 
and that this result is independent of the SED used. Using 
the calibration by Vacca et al. (1996), this corresponds 
to a spectral type earlier than 07.5. This calibration is. 



however, based on plane-parallel models which do not in- 
corporate stellar winds and line blocking/blanketing. A 
more recent calibration by Martins et al. (2002), based 
on non-LTE line blanketed atmosphere models with stel- 
lar winds computed with the CMFGEN code of Hillier & 
Miller (1998), yields a spectral type equal to or earlier 
than 06. This spectral type determination is compared to 
estimates using other methods in Sect. 5. 



5. The ionizing star of G29.96 

The opening up of the infrared window has favoured in re- 
cent years the use of infrared spectra to infer the spectral 
type of the ionizing stars in heavily extincted regions. The 
first spectral classification in the K-band of a star ionizing 
an UCHii region was that of G29.96 (Watson & Hanson, 
1997). The K-band classification (Hanson et al., 1996) is 
based on the presence and strength of photospheric lines of 
He I (2.058 and 2.1126 /xm),Civ (2.078 /im),Niii (2.1155 
(um). Hi (2.1661 ^im) and Hen (2.1885 ^m). The Hi 
and He i lines, however, arc heavily contaminated by neb- 
ular emission. Based on the other features, the spectral 
types that can be determined from the direct observation 
of the K-band spectrum are: 03-04 (N ill emission and 
Hen absorption, but no Civ emission), 05 06.5 (both 
Niii and C IV emission and Hell absorption), 07 08 
(N III emission and He ii absorption, and weak C iv emis- 
sion), and 09 or later (none of these). Using this classi- 
fication scheme, Watson & Hanson (1997) restricted the 
ionizing star in G29.96 to spectral type 05-06.5. 

Figure 11 shows the spectrum of the ionizing star 
of G29.96 extracted from the slit (Kaper et al. 2002*; 
Bik et al. in prep.). The emission features of CiV and 
Niii are clearly present. The spectral coverage of our ob- 
servations was not large enough to include the He ii line 
at 2.1885 /im. The presence and strength of the C iv and 
N III features is consistent with the estimate by Watson 
& Hanson (1997), i.e. with an 05 06.5 star. Watson & 
Hanson (1997) note that spectral types as late as 07-08 
could also be consistent with the data if the C IV and 
Niii features were enhanced as a consequence of G29.96 
being a region with a metallicity twice that of the Sun. 
However, these late spectral types are ruled out because 
they are not hot enough to produce a He+ zone coinci- 
dent with that of H+ (see Sect. 4.2). Thus, the K-band 
spectrum of the ionizing star, together with the condition 
imposed by the extent of the He"*" zone (equal to or earlier 
than 06), limits the spectral type to 05-06. 

Various other estimates of the spectral type of the 
dominant ionizing star of G29.96 have been obtained from 
very different observations and methods. We briefly dis- 
cuss these estimates using up-to-date stellar models and 
homogeneous assumptions. In particular, all observational 



We note that Kaper et al. wrongly identify the He: 
2.1649 //m line as Hen 2.1652 yum, and use this identificar 
tion to argue that the ionizing star is of even earlier spectral 
type than 05. 
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Method 


Teff (kK) 


Spectral type 


Reference 


Bolometric luminosity 


< 


< 03* 


1, 2 


Lyman continuum photon flux 


35 - 40t 


08 - 06* 


3, 4, 5, 6, 7 


Ncar-IR photometry 


< 42.5 


< 05* 


2 


Nebular IR lines 


32 - 38 


09 - 06.5* 


8 


Extent of the He"*" zone 


> 40" 


> 061' 


9 


K-l)aii(l spectrum 


42..")* 


0G.5 0.^) 


9.10 



( ) ">" means "earlier than", and "<" means "later than". (^) Adopting the calibration by Vacca ct al. (1996). (*) Adopting 
the Toff scale by Martins et al. (2002) for dwarf stars. (") See Sect. 4.2. References: (1) Peeters et al. (2002); (2) Watson et al. 
(1997); (3) Wood & Churchwell (1989); (4) Cesaroni et al. (1994); (5) Becker et al. (1994); (6) Fey et al. (1995); (7) Kim & Koo 
(2001); (8) Morisset et al. (2002); (9) this work; (10) Watson & Hanson (1997); 
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Fig. 11. Spectrum of the ionizing star in G29.96. The emis- 
sion features of Civ at 2.080 /itm and Niii at 2.116 /xm are 
indicated. These features constitute the base of the K band 
classification. The lines of Hi and He I are contaminated by 
nebular emission. 



constraints are scaled to a distance of 6 kpc and we use 

the latest temperature scale of O stars by Martins ct al. 
(2002). We assume that the ionizing star is an unevolved 
zero age main sequence (ZAMS) star, although we note 
that some authors (Watson ct al., 1997; Morisset ct al., 
2002) suggest that it may be a more evolved star. These 
estimates are (see also Table 2): 

• Bolometric luminosity The total bolometric luminos- 
ity obtained from the IRAS fluxes (Peeters et al., 
2002) and the overall SED (Watson et al., 1997) in 
an arcminute-sized region is log{L/LQ) = 5.9. Since 
this luminosity is integrated over a large region, it is 
likely to inc;lu(lc contributions from other sources, for 
instance, the hot core. Hence, it must be considered 
as an upper limit to the bolometric luminosity of the 
ionizing star. This upper limit to the luminosity cor- 
responds to a star with an effective temperature lower 
than 49000 K (Vacca et al., 1996). Adopting the T^s- 



spcctral type calibration by Martins ct al. (2002) for 
dwarf stars, this corresponds to a star later than 03. 

• Number of Lyman continuum photons The number of 
Lyman continuum photons derived for G29.96 from ra- 
dio continuum observations is in the range log(AfLyc) = 
48.43-48.93 s'^ (Wood & Churchwell, 1989; Cesaroni 
ct al., 1994; Becker et al., 1994; Fey ct al., 1995; Kim 
& Koo, 2001), which corresponds to an effective tem- 
perature in between 35 000 and 40000 K (Vacca et al., 
1996). This Tcs range corresponds to a star with a 
spectral type of 08-06 (Martins et al., 2002). 

• Near-infrared photometry Constraints on the K mag- 
nitude and the bolometric luminosity allowed Watson 
et al. (1997) to place limits on the location of the 
ionizing star in the H-R diagram. They found a 3a 
upper limit on the effective temperature of 42 500 K. 
Adopting the Teg-spectral type calibration by Martins 
et al. (2002), this corresponds to a star equal to or 
later than 05. 

• Nebular infrared lines The ratios of nebular fine- 
structure lines ([Ariii]/[Arii] 9.0/7.0, [Siv]/[Siii] 
10.5/18.7, [Nciii]/[Ncn]15. 5/12.8 and [Niii]/[Nii] 
57/128 A<m) observed by ISO (Peeters et al., 2002) 
have been used by Morisset et al. (2002) to constraint 
the properties of the ionizing star in G29.96. They use 
the most recent non-LTE line blanketed atmospheres 
with stellar winds and obtain effective temperatures in 
the range 32 000-38000 K. Adopting the Teff-spectral 
type calibration by Martins et al. (2002), this corre- 
sponds to a star with a spectral type 06.5-09. 

Practically all the indirect methods to derive the 
spectral type which have been listed above agree well 
with the spectral type derived from the K-band spec- 
trum. The exception is the method based on the nebu- 
lar infrared lines, which predicts a slightly later type. In 
Martm-Hernandez et al. (2002b) it was indicated that the 
infrared fine-structure line ratios such as [Neiii]/[Neii] 
15.5/12.8 /xm are influenced by the stellar and nebular 
metallicity. Hence, their diagnostic use as indicator of the 
degree of ionization of the nebula is conditioned to an 
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adequate treatment of the effect of metalUcity on the stel- 
lar UV spectrum, i.e, on the line blocking/blanketing and 
stellar wind. 

6. The dynamics of G29.96 

Any model for the formation and evolution of the UCHii 
region G29.96 should be able to explain the following char- 
acteristics: 

1. The general morphology of the nebula is described by a 
bright arc with the ionizing star in its focus and a lower 
surface brightness tail (see the 2 cm radio continuum 
image of Wood & Churchwell 1991, as well as the K- 
band image of Lumsden & Hoarc 1996). 

2. The molecular gas traced by H2 shows a similar bright 
arc (see Fig. 1), but with the peak displaced 'deeper' 
into the molecular cloud (by 0'.'93) than the ionized 
gas. 

3. The ionized gas shows a steep gradient (see Fig. 4) 
which goes from —31 km (just in front of the bright 
arc) to +7 kms~^ (just behind the position of the 
star) relative to the molecular cloud velocity (which is 
about 98 kms~^). The velocity then drops again to 
—4 kms"""^ in the tail. 

4. In contrast with the large variations in velocity ob- 
served for Br7 across the symmetry axis of G29.96, 
the H2 emission shows no significant velocity gradients. 
Indeed, it seems to be at the velocity of the molecular 
cloud. 

5. The FWHM of the Br7 and H2 lines are significantly 
higher in front of the respective bright arcs than in the 
tail. Specifically, the linewidths observed for Br7 range 
from 42 — 46 kms~^ behind the star to values as high 
as 62 km s~^ just in front of the arc. In the case of the 
H2 line, the linewidth increases from 33 kms~^ in the 
tail to ~ 50 kms^^ in front of the H2 arc. 

Two different models aimed at explaining the forma- 
tion and evolution of cometary H II regions and applied to 
G29.9 appear in the literature: the bow shock model and 
the champagne flow model. Neither of these two models 
can account for all the characteristics observed for G29.96. 
We critically discuss both models. 

6.1. The stellar bow shock model 

The morphology of cometary Hll regions, in particular 
their near-parabolic shape, has been interpreted in terms 
of a bow shock (Wood & Churchwell, 1989). In such a 
model, an OB star moves with supersonic velocity (of a 
few kms~^) inside a molecular cloud. The interaction of 
the stellar wind with the ambient medium results in the 
formation of a bow shock with a cometary morphology. 
Detailed analytical and numerical models have been devel- 
oped and applied to the UCHii region G29.96 (van Buren 
et al., 1990; Mac Low et al., 1991; van Buren & Mac Low, 
1992). 



The structure proposed by van Buren et al. (1990) of 
an UCHii region trapped in a bow shock is the following. 
The fast stellar wind (with a terminal velocity in excess 
of 2000 kms~^) rams into the molecular cloud material, 
which is moving at supersonic velocity with respect to the 
star, resulting in the formation of a shock wave. The shock 
occurs at a radius where the momentum flux in the wind 
equals the ram pressure of the ambient medium. The im- 
mediate post-shock gas temperature is ~ 10^ K. At this 
temperature, radiative cooling is not very efficient, but 
cooling via thermal conduction can lower the temperature 
very rapidly and consequently, this layer of shocked wind 
gas is rather thin (Comeron & Kaper, 1998). As the star 
moves through the ambient molecular medium, the hot 
bubble of shocked wind gas sweeps up material into a thin 
parabolic shell. Because of the rapid cooling, the wind only 
transfers its momentum to the swept-up shell in what is 
called a "momentum-driven" bow-shock. This shell traps 
all or part of the ionizing flux of the star, forming a thin 
layer of dense ionized gas surrounded by shocked molec- 
ular gas. The layer of ionized gas has the properties of a 
typical radiatively excited H 11 region at a temperature of 
~ 10^ K. 

In terms of the dynamics, the bow shock models essen- 
tially predict a redshift or blueshift of the ionized gas rela- 
tive to the ambient cloud material depending on whether 
the star is moving away from the observer and into the 
molecular cloud, or towards the observer. The material 
in the bow shock is transported (due to a pressure gra- 
dient) from the head to the tail with a velocity that will 
not exceed the space velocity of the star by much. Note, 
however, that the material inside the bow shock flows in 
a direction away from the head to the tail, i.e. at some 
point it moves in a direction opposite from the stellar mo- 
tion. The maximum velocity difference is then of the order 
of Vi,cos{i), where v^, is the relative velocity of the star 
through the ISM and i is the inclination angle. As empha- 
sized by Lumsden & Hoare (1996), the velocity structure 
observed for the Br7 line can only be explained by the bow 
shock model if the molecular cloud velocity is 67 kms^^, 
rather than the 98 kms^^ obtained from rotational lines 
(e.g. Churchwell et al., 1990, 1992; Cesaroni et al., 1992; 
Bronfman et al., 1996: Olmi & Cesaroni, 1999). Moreover, 
if we attribute the difference of ^ 32 kms~^ observed in 
between positions — 7"and 0" (see Fig. 4 a) to the star- 
cloud interaction, then, with the inclination angle esti- 
mated from the morphology of —45 degrees (van Buren 
& Mac Low, 1992), this translates into a stellar velocity 
of ~ 45 kms~^. This is a very high velocity for a star 
in a cluster. Now, the difficulty with the stellar velocity 
might be mitigated if the velocity shoulder from 67 to 87 
kms~^ observed in Fig. 4 a between positions —7" and 
—3" traces motions in the hot core located just 2" west of 
the bright arc. This hot core is known to have an embed- 
ded protostellar object with a prominent outflow (A. Gibb, 
private communication). In that case, the required molec- 
ular cloud velocity is 87 kms~^ and the stellar velocity 
needed to explain the velocity difference between positions 
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— 3"and 0" is only 17 kms"^. However, still this velocity 
of 87 kms"^ is much less than the observed molecular 
cloud velocity of 98 kms"-'^. 

The bow shock model provides a natural explanation 
for the general morphology of cometary H II regions if the 
inclination angle is larger than ±30 degrees (Mac Low 
et al., 1991; van Buren & Mac Low, 1992). As an aside, 
we can test whether the morphology observed for G29.96 
is consistent with that of an 05-06 main sequence star 
(see Sect. 5) developing a bow shock. Van Buren et al. 
(1990) propose a simple analytical model for bow shocks 
in molecular clouds. The standoff distance Isd, i.e. the dis- 
tance from the star where the momentum flux in the wind 
equals the ram pressure of the ambient medium, is given 
by: 



«sd(pc) = 0.015 



lO-6M0/yr 



1/2 



10^ cm/s 



1/2 



/ riH / 
Vl05cm-3; V 



lO^cm/s 



where rh* is the stellar wind mass-loss rate, Ww is the ter- 
minal velocity of the wind, nn is the number density of 
hydrogen nuclei in all forms in the ambient gas and is 
the relative velocity of the star through the ISM. Wo have 
adopted a mean mass per particle, /in. of 1.4 for the molec- 
ular gas. We take an 06 star with a inetallicity twice solar 
(see Sect. 4.2) and stellar parameters adopted from Vacca 
et al. (1996). Following Vink et al. (2001), the correspond- 
ing stellar wind parameters are m^, = 2.81 x 10~^ Mq yr~^ 
and Ww = 2951 kms~^. A star of these characteristics 
moving supersonically trough the molecular cloud at a ve- 
locity of 17 kms^^ will develop a bow shock with a stand- 
off distance of 0.025 pc for a typical ambient gas density 
of 10^ cm~^. At a distance of 6 kpc and for an inclina- 
tion angle of —45 degrees, this distance corresponds to a 
projected separation of 0'.'6. From the observations (see 
Sect. 3), the projected standoff distance would be the off- 
set between the star and the peak of the Br7 arc, i.e. 2'.'27. 
The observations imply either a much lower stellar veloc- 
ity (of the order of 5 kms~^) or a molecular gas density 
at the stagnation point below 10* cm~^. 

6.2. The champagne flow model 

G29.96 has also been interpreted in terms of a champagne 
flow. Champagne flow (blister) models assume that the 
medium in which a massive star is born is not uniform but 
has strong density gradients, which can give rise to an H II 
region that expands supersonically away from the high- 
density region in a so-called champagne flow (e.g. Tenorio- 
Tagle, 1979; Bodenheimer et al., 1979). Theoretical simu- 
lations of the expansion of an H II region in the champagne 
phase, i.e. when the ionization front of the nebula crosses 
a region of strong density gradients, suc;li as the edge of 
the molecular cloud, and expands into the lower density 
intercloud medium, have been presented by Yorke et al. 



(1983). During this champagne flow phase, the Hii region 
is ionization bounded on the high-density side (i.e. towards 
the molecular cloud) and density bounded on the side of 
the outward champagne flow. 

The morphology of the ionized gas as predicted by 
these models is characterized by a bright compact compo- 
nent and an extended low brightness component roughly 
shaped like an opened fan. This compact component is 
inconsistent with the morphology observed for G29.96. At 
the same time, the opening of the morphology predicted 
for the tail region is not evident in G29.96. 

In terms of the velocity structure, the numerical cal- 
culations (Yorke et al., 1983) show that the ionized gas is 
accelerated in the direction away from the molecular cloud 
to a relatively high velocity, which can attain values in ex- 
cess of 30 kms~^. On the contrary, very little variations 
should be observed near the core of the H ii region, where 
the ionization front should be essentially at the molecular 
cloud velocity (98 kms~^ in the case of G29.96). However, 
we see it at 67 kms~^ (or 87 kms~^ if we attribute the ve- 
locity shoulder between position —7" and —3" to velocity 
flows within the hot core, see Sect. 6.1). 

6.3. Final remarks 

Neither of the two competing models can explain the ob- 
served velocity structure for G29.96. Partly, this may re- 
flect the complex structure of this region, where the pres- 
ence of a hot core with an embedded young stellar object 
and an associated outflow may obfuscate the dynamics. 

Further progress will require the direct determination 
of the stellar velocity from photospheric absorption lines, 
for instance, from near-IR spectroscopy at high spectral 
resolution. We tried to obtain information on the stellar 
velocity from the C IV and N III emission lines observed in 
our K-band spectrum, but it was not possible due to the 
complexity of these lines together with the limited spectral 
resolution (Av ~ 33 kms~^). Also, a complete mapping 
of the dynamics of the ionized and molecular gas at both 
high spectral and spatial resolution will be essential. 

A comparison of the observations with alternative 
models will also be important. Several authors, for in- 
stance, have investigated the interaction of a stellar wind 
and a clumpy molecular cloud material (e.g. Dyson et al., 
1995; Redman et al., 1996; Williams et al., 1996). In this 
mass-loaded model, material from dense neutral clumps of 
gas is added to the ionized region around a newly formed 
massive stars; this neutral material soaks up ionizing pho- 
tons, slowing down the expansion of the UCHii region. 
This mass-load model can explain cometary regions if the 
star is located in a density gradient of mass-loading clumps 
(Redman et al., 1998). G29.96 provides an excellent labo- 
ratory to test this model. 

7. Summary 

We presented a high quality, medium resolution K-band 
spectrum of the ultracompact region G29.96— 0.02 ob- 
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tained with the VLT. The sht was positioned along the 
symmetry axis of the nebula. Besides the spectrum of 
the embedded ionizing star, the observations included the 
emission-line spectrum produced by the ionized gas with 
sub-arcsec spatial resolution. The nebular spectrum in- 
cludes Br7, several helium emission lines, and a molecular 
hydrogen line. 

The study of the relative strength of the Br7 and He i 
lines as a function of position along the slit allowed us 
different lines of investigations. First, the observed he- 
lium lines in the nebula have been confronted with the 
He I recombination model. We have shown that the de- 
viations of the transition of 4'^S-3'^P at 2.1128 /xm from 
the Case B approximation are well described by the cur- 
rent understanding of the collisional and radiative trans- 
fer effects on the He i recombination spectrum of nebu- 
lae. The Case B approximation reproduces the ratios of 
He I lines originating from higher quantum states reason- 
ably well. These Hei line ratios indicate a constant elec- 
tron temperature, between 6400 and 7500 K. Second, the 
relative strength of the Br7 and the singlet transition of 
4^S-3ip at 2.1137 /xm indicates that He+/H+ is practi- 
cally constant along the slit, with an average abundance of 
0.087 ± 0.001; we argue that He is singly ionized through- 
out the nebula. This conditions sets a lower limit to the 
effective temperature of the star. 

The K-band spectrum of the ionizing star, together 
with the constraint imposed by the coincidence of the He+ 
and H+ zones, limit its spectral type to 05-06. The vari- 
ous observational constraints on the effective temperature 
of the ionizing star are reviewed, and they are shown to 
be in rather good agreement with the above spectral type. 

Finally, the variations in peak velocity and width of 
the Br7 and H2 lines allowed us to study the kinematical 
structure of the ionized and molecular gas at higher spatial 
resolution than in previous studies. We find that neither 
the wind bow shock model nor the champagne flow model 
are supported by the observations. A firm constraint on 
the velocity of the star and a complete mapping of the 
dynamics of the ionized and molecular gas at both high 
spectral and spatial resolution will be essential to under- 
stand the formation and evolution of G29.96. 

Acknowledgements. We would like to thank the referee. Dr. P. 
Conti. MICE is supported at MPE by DLR (DARA) under 
grants 50 QI 86108 and 50 QI 94023. 

References 

Afflerbach, A., Churchwell, E., Acord, J. M., et al. 1996, 

ApJS, 106, 423 
Afflerbach, A., Churchwell, E., Hofner, P., & Kurtz, S. 

1994, ApJ, 437, 697 
Becker, R. H., White, R. L., Hclfand, D. J., & 

Zoonematkermani, S. 1994, ApJS, 91, 347 
Benjamin, R. A., Skillman, E. D., & Smits, D. P. 1999, 

ApJ, 514, 307 
— . 2002, ApJ, 569, 288 



Bodenheimer, P., Tenorio-Tagle, G., & Yorke, H. W. 1979, 
ApJ, 233, 85 

Bragg, S. L., Smith, W. H., & Brault, J. W. 1982, ApJ, 

263, 999 

Bronfman, L., Nyman, L. ., & May, J. 1996, A&AS, 115, 

81 

Cesaroni, R., Churchwell, E., Hofner, P., Walmsley, C. M., 

& Kurtz, S. 1994, A&A, 288, 903 
Cesaroni, R., Hofner, P., Walmsley, C. M., & Churchwell, 

E. 1998, A&A, 331, 709 
Cesaroni, R., Walmsley, C. M., & Churchwell, E. 1992, 

A&A, 256, 618 
ChurchweU, E. 1990, A&AR, 2, 79 

Churchwell, E., Walmsley, C. M., & Cesaroni, R. 1990, 

A&AS, 83, 119 
ChurchweU, E., Walmsley, C. M., & Wood, D. O. S. 1992, 

A&A, 253, 541 
Comcron, F. & Kapcr, L. 1998, A&A, 338, 273 
De Buizer, J. M., Watson, A. M., Radomski, J. T., Piiia, 

R. K., & Telesco, C. M. 2002, ApJ Lett., 564, LlOl 
Deharveng, L., Peha, M., Caplan, J., & Costero, R. 2000, 

MNRAS, 311, 329 
Dyson, J. E., Williams, R. J. R., & Redman, M. P. 1995, 

MNRAS, 277, 700 
Ferland, G. J. 2000, Hazy, a brief introduction to 

Cloudy 94.00a (Department of Physics and Astronomy, 

University of Kentucky) 
Fey, A. L., Gaume, R. A., Claussen, M. J., & Vrba, F. J. 

1995, ApJ, 453, 308 
Garay, G. & Lizano, S. 1999, PASP, 111, 1049 
Hanson, M. M., Conti, P. S., & Riekc, M. J. 1996, ApJS, 

107, 281 

Hanson, M. M., Luhman, K. L., & Rieke, G. H. 2002, 

ApJS, 138, 35 
Hillier, D. J. & Miller, D. L. 1998, ApJ, 496, 407 
Hofner, P. & Churchwell, E. 1996, A&AS, 120, 283 
Hummer, D. G. & Storey, P. J. 1987, MNRAS, 224, 801 
Kaper, L., Bik, A., Hanson, M., & Comcron, F. 2002, in 

Hot Star Workshop HI: The Earliest Stages of Massive 

Star Birth; ed. P.A. Crowther, ASP Conference Series, 

267, 95 

Kim, K. & Koo, B. 2001, ApJ, 549, 979 
Kingdon, J. & Ferland, G. J. 1995, ApJ, 442, 714 
Lumsden, S. L. & Hoare, M. G. 1996, ApJ, 464, 272 
— . 1999, MNRAS, 305, 701 

Lumsden, S. L. & Puxley, P. J. 1996, MNRAS, 281, 493 
Lumsden, S. L., Puxley, P. J., & Hoare, M. G. 2001a, 

MNRAS, 320, 83 
— . 2001b, MNRAS, 328, 419 

Mac Low, M., van Buren, D., Wood, D. O. S., & 

Churchwell, E. 1991, ApJ, 369, 395 
Martm-Hernandez, N. L., Peeters, E., Morisset, C, et al. 

2002a, A&A, 381, 606 
Martm-Hernandez, N. L., Vermeij, R., Tielens, van der 

Hulst, J. M., & Peeters, E. 2002b, A&A, 389, 286 
Martins, F., Schaerer, D., & Hiller, D. J. 2002, A&A, 382, 

999 



N.L. Martm-Hernandez et al.: A VLT spectroscopic study of the ultracompact Hii region G29.96— 0.02 



15 



Maxia, C, Testi, L., Cesaroni, R., & Walmsley, C. M. 

2001, A&A, 371, 287 

Morisset, C, Schaerer, D., Martm-Hernandez, N. L., et al. 

2002, A&A, 386, 558 

Olmi, L. & Cesaroni, R. 1999, A&A, 352, 266 
Pauldrach, A. W. A., Hoffmann, T. L., & Lennon, M. 

2001, A&A, 375, 161 

Peeters, E., Martm-Hernandez, N. L., Damour, F., et al. 

2002, A&A, 381, 571 

Pratap, P., Megeath, S. T., & Bergin, E. A. 1999, ApJ, 
517, 799 

Pratap, P., Menten, K. M., & Snyder, L. E. 1994, ApJ 

Lett., 430, L129 
Redman, M. P., Williams, R. J. R., & Dyson, J. E. 1996, 

MNRAS, 280, 661 
— . 1998, MNRAS, 298, 33 

Schaerer, D. & de Koter, A. 1997, A&A, 322, 598 
Shaver, P. A., McGee, R. X., Newton, L. M., Banks, A. C, 

& Pottasch, S. R. 1983, MNRAS, 204, 53 
Tenorio-Tagle, G. 1979, A&A, 71, 59 
Vacca, W. D., Garmany, C. D., & ShuU, J. M. 1996, ApJ, 

460, 914 

van Buren, D. & Mac Low, M. 1992, ApJ, 394, 534 

van Buren, D., Mac Low, M., Wood, D. O. S., & 

ChurchweU, E. 1990, ApJ, 353, 570 
Vanzi, L., Rieke, G. H., Martin, C. L., & Shields, J. C. 

1996, ApJ, 466, 150 
Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, 

A&A, 369, 574 
Walsh, A. J., Burton, M. G., Hyland, A. R., & Robinson, 

G. 1998, MNRAS, 301, 640 
Watson, A. M., Coil, A. L., Shepherd, D. S., Hofner, P., 

& ChurchweU, E. 1997, ApJ, 487, 818 
Watson, A. M. & Hanson, M. M. 1997, ApJ Lett., 490, 

L165 

WiUiams, R. J. R., Dyson, J. E., & Redman, M. P. 1996, 

MNRAS, 280, 667 
Wood, D. O. S. & ChurchweU, E. 1989, ApJS, 69, 831 

— . 1991, ApJ, 372, 199 

Yorke, H. W. 1986, ARA&A, 24, 49 

Yorke, H. W., Tenorio-Tagle, G., & Bodenheimer, P. 1983, 
A&A, 127, 313 



